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We study heat conduction in one dimensional lattice dynamical systems far from equilib- 
rium. The Fermi-Pasta-Ulam model and the tfi 4 model are numerically compared to elucidate 
differences between momentum-conserving and nonconserving systems. As a results, it is 
found that the heat flux in the </> 4 model does not increase monotonically as the temperature 
differences at the ends of the lattice is increased, while it does in the FPU chain. 
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Currently, heat conduction in one-dimensional systems is getting more and more attention 
in nanoscale physics. Murayama studied heat conduction in the single walled carbon nanotube 
and found that the thermal conductivity behaves like one-dimensional systems. 1 Namely, the 
thermal conductivity a nanotube with a small diameter did not converge to a finite value with 
increase in tube length, but obeyed a power law relation as is observed in the one-dimensional 
systems. Similar behavior was observed in the nanoscale world by several authers. 2,3 

In macroscopic systems, heat conduction obeys the Fourier's law 

J = -kVT, (1) 

where J and VT denote the heat flux and the temperature gradient. Recent studies in one- 
dimensinal systems have clarified that size-dependence of the thermal conductivity n is dif- 
ferent between momentum-conserving and nonconserving systems. For example, the thermal 
conductivity diverges as k ~ iV ' 37 in the FPU chain, where the momentum is a conserved 
quantity. 4 On the other hand, it converges in the limit of large system size in the <f> model, 
where momentum is not conserved due to the existence of on-site potentials. This principle 
is also applied to the diatomic Toda chain, 11 the Frenkel-Kontrova model, 6 ' the ding-a-ling 
model, 8 the ding-dong model, 9 and many others, though the integrable systems like the har- 
monic chain or the Toda lattice show another behavior k ~ N and there are some exceptions 
to the principle. 12,14 Now, this behavior is understood as an effect of long-time tail in the 
autocorrelation function of the total heat flux. 13 
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On the other hand, in nanoscale devices, the temperature gradient need not be small 
and heat conduction is not necessarily described by the Fouries's law. Thus, heat conduction 
beyond the linear-response regime can be observed in nanoscale physics. Accordingly, we study 
in this paper properties of heat conduction far from equilibrium, and particularly look into 
differences between momentum-conserving and nonconserving models. 

We consider one-dimensional systems with Hamiltonian of the type 

N ,- 2 1 N 

H = E y + U M + E v ^ ~ *). (2) 



i=i 



i=0 



where qi and pi denote the displacement and the momentum of the ith particle, U(qi) the 
on-site potential and V(qi+i —qi) the interaction potential between nearest-neighbor particles. 
We use the fixed boundary condition, whch is represented by putting qo = qN+i = in Eq. 
(2). We consider the following two models: (i) the FPU model is described by V(q) = \q 2 + \q A 



and U{q) = and (ii) the </> 4 model is decribed by V(q) = \q 2 and U(q) = \q 2 + \q A . 



2^ 1 4 L - 

Moreover, two heat reservoirs are attached to the system, that is, particle 1 is in contact 
with a heat reservoir at temperature Tl and particle N is in contact with another heat 
reservoir at temperature Tr. In the present study, we use three kinds of heat reservoirs, the 
Langevin thermostats, the Nose-Hoover thermostats, and the thermal wall. In the Langevin 
thermostats, random forces and dissipation terms are introduced into the equations of motion. 
Namely, the equations of motion of particles 1 and ./V are 

Pi = -U'(qi) + V'(q 2 - Ql ) - V'( qi ) - 7LP1 + 

Pn = -U'(q N ) + V'(-q N ) - V'{q N - q N -i) - IrPn + £,n(t), 

where j a denotes the coupling strength between the particle and the heat bath and £l (t) and 
£r(£) are the Gaussian white noise with zero mean and variance 

(UtW)) = 2 la T a S a( sS(t - 0, (a,P = R,L). (3) 

In the case of Nose-Hoover thermostats, deterministic forces play the role of heat reservoir. 
Namely, the equations of motion for particles 1 and ./V are modified as 

Pi = -U'( qi ) + V'(q 2 - Ql ) - V'( qi ) - to, ( L = 1^-1^, 

PN = -U'(q N ) + V'{-q N ) - V'(q N - q N ^) - ^PN, U = ^ ~ lj , 

where is the thermostat response time. In the thermal wall, if a particle hits the wall, it is 
reflected with a momentum randomly chosen according to the distribution 

/M = Mexp(-0 (4) 
where T is the temperature of the wall. 
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Fig. 1. Temperature profile for the FPU model (a) and the </> 4 model (b) with system size N=32 and 
the Langevein thermostats, the Nose-Hoover thermostats and the thermal walls at Tl — 2.0 and 
T R = 1.0. 



In our numerical simulations, the 4th-order Runge-Kutta and the 6th-order symplectic 
integrator were used in the Langevin and the Nose-Hoover cases. In the cases of the thermal 
walls, the 4th-order Symplectic Integrator was employed. All the simulations were done with 
a time step At = 0.01 and the parameters are chosen as jl = 1R = 0.5 and 6 = 0.1. 
The temperature of particle i is defined as kinetic temperature Tj = (pf) and local heat flux 
through particle i is given by 

3i = \{Pi+i +Pi) F (Qi+i ~ Qi)- (5) 

Both the quantities were obtained with averaging over 10 10 time steps in the 4> 4 model and 
10 9 time steps in the FPU model after a nonequilibrium steady state was reached. 

Figure 1 shows the temperature profile for the 4 model and the FPU model near equilib- 
rium steady states. We see that smooth temperature profile is formed except for boundaries, 
where temperature gaps exist. 

When the temperature difference between the two heat reservoirs is larger, the system is 
beyond the linear response regime and then Fourier's law does not hold. Figure 2 shows the 
temperature profile in the model and the FPU model, respectively. In the FPU model, 
the temperature profile is qualitatively similar to the near-equilibrium case. Only difference 
is that the gaps between the heat reservoirs and particle 1 or N become larger than the near- 
equilibrium case. On the other hand, in the 4> 4 model, the temperature profile is drastically 
transformed. Namely, local temperatures are largely shifted to low values in the bulk region 
and changes are notable only near the end at the higher temperature. 

Heat flux also shows peculiar behavior. Figures 3 and 4 show the mean heat flux as a 
function of the temperature difference Xjv — T\. 

1 N 

i=i 
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Fig. 2. Temperature profile for the FPU model (a) and the <fi 4 model (b) with system size N=32 and 
the Langevein thermostats, the Nose-Hoover thermostats and the thermal walls at Tl = 40.0 and 
T R = 1.0. 





Fig. 3. Averaged heat flux and temperature difference at both ends for the FPU model with the 
Langevein thermostats (a), the Nose- Hoover thermostats (b) and the thermal walls (c). 



where ji is the local heat flux at the ith site defined by 

ji = \{Pi+l +Pi) F (Qi+l ~ Qi), (7) 

where F(x) = — ^V(x). It is observed that as the temperature difference increases, the 
averaged heat flux increases monotonically in the FPU model. In the </> 4 model, however, the 
averaged heat flux saturates at some value of temperature difference and does not increase 
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Fig. 4. Averaged heat flux and temperature difference at both ends for the </> 4 model with the 
Langevein thermostats (a) , the Nose- Hoover thermostats (b) and the thermal walls (c) . 



further. Similar results are obtained irrespective of the types of heat reservoirs. Therefore, 
these differences are caused only by the characteristic of the lattice dynamical systems. 

Moreover, we compute the harmonic interaction potential with and without harmonic 
on-site potential cases. In these cases, we confirmed that the averaged heat flux increases 
monotonically. Thus, the peculiar behavior of the (j) 4 model should be attributed to com- 
bination of nonintegrability and lack of momentum conservation. Extensive study to clarify 
the origin of our findings is planned for the future. Our results may provide a guideline for 
designing nanoscale devices. 

We have studied the heat conduction far from equilibrium in non equilibrium steady states 
in the FPU and (p 4 chains. As a results, we have found the new interesting phenomenon of 
the heat conduction. In the FPU model, heat flux increases monotonically as temperature 
difference between particles 1 and N. However, in the <p 4 model, heat flux does not increase 
monotonically. This property does not depend on the heat reservoirs. 
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